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Abstract 

This talk presents the results of our study of systematic errors of the ground- 
state parameters obtained by Shifman-Vainshtein-Zakharov (SVZ) sum rules. 
We use the harmonic-oscillator potential model as an example: in this case we 
know the exact solution for the polarization operator, which allows us to obtain 
both the OPE to any order and the parameters (masses and decay constants) of 
the bound states. We extract the parameters of the ground state by making use 
of the standard procedures of the method of QCD sum rules, and compare the 
obtained results with their known exact values. We show that if the continuum 
contribution to the polarization operator is not known and is modelled by some 
effective continuum threshold, the standard procedures adopted in sum rules 
do not allow one to gain control over the systematic errors of the extracted 
ground-state parameters. 



A QCD sum-rule calculation of hadron parameters involves two steps: 
one first constructs the operator product expansion (OPE) series for a relevant 
correlator and then extracts the parameters of the ground state by a numerical 
procedure. Each of these steps leads to certain uncertainties in the final result. 

The first step lies entirely within QCD and, in the case of SVZ sum rules, 
allows for a rigorous treatment of the uncertainties: the correlator is not known 
precisely because of uncertainties in quark masses, condensates, as, etc., but all 
corresponding errors in the correlator may be controlled. [Complications arising 
in light-cone sum rules are discussed in our second talk 2).] 

The second step lies beyond QCD: even if several terms of the OPE for the 
correlator were known precisely, the hadronic parameters might be extracted by 
a sum rule only within some error, which may be treated as a systematic error of 
the method. 

Here we present the results of our recent study of systematic uncertainties 
of the sum-rule procedures!^ , Xo this end, a quantum-mechanical harmonic- 
oscillator (HO) potential model is a perfect tool: in this case both the spectrum 
of bound states (i.e., masses and wave functions) and the exact correlator (and 
hence its OPE to any order) are known precisely. Therefore, one may apply the 
sum-rule machinery for extracting parameters of the ground state and test the 
accuracy of the extracted values by comparing with the known exact results. In 
this way the accuracy of the method can be probed. For a detailed discussion of 
various aspects of sum rules in quantum mechanics, we refer to Refs. [5-9]. 

To illustrate the essential features of the QCD calculation, we consider a 
non-relativistic model with a confining potential, 
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y(r) = ^— , r=|r|, (1) 

and analyze the Borel transform H(/i) of the polarization operator H(£'), which 
gives the evolution operator in the imaginary time l/fj.: 
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For the HO potential ((T]), the exact analytic expression for n(pt) is well known: 



Expanding the above expression in inverse powers of fi, we get the OPE series 
Hope (m) = (a^) + Hi (^) + Ha (m) + • ■ ■ 



higher power corrections may be derived from the exact result 

The "phenomenological" representation for n(/i) is obtained by using the 
basis of hadronic eigenstates of the model, namely, 

°° / E \ 
n(M) = ^i?„exp ^ , (5) 
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where En is the energy of the nth bound state and i?„ [the square of the 
leptonic decay constant of the nth bound state] is given by 
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Ra=[ — ] |*«(r - 0)1^. (6) 



For the lowest states, one finds from ^ 

Eo = ^w, Ro ^ 2V2u;^^^, Ei = ^lu, Ri ^ 3\/2cj3/^ .... (7) 

The sum rule is just the equality of the correlator calculated in the "quark" basis 
and in the "hadron" basis: 
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Following Ref. [1], we use explicit expressions for the power corrections, but for 
the zeroth-order free-particle term we use its expression in terms of the spectral 
integral. 

Let us introduce an "effective" continuum threshold Zcsip), different from 
the physical /i- independent continuum threshold Zcont, by the relation 
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ncont(Ai) = J dz pphcn{z) exp ^ J dz po{z) exp - 

Zcont ZctK/j) 



(9) 



The spectral densities Pphcn(-z) and po{z) are different functions. Thus the two 
sides of © can be equal to each other only if the effective continuum threshold, 
ZcsifJ-), depends on /i in an appropriate way. In our model, we can calculate 
Ilcont precisely and, therefore, we can obtain the function Zoff (m) by solving ([5]). 
In the general case of an actual QCD sum-rule analysis, the effective continuum 
threshold is not known and constitutes one of the essential fitting parameters. 
Making use of we now rewrite the sum rule ^ in the form 

iJoexp ( ) = n(^, Zoff(Ai)), (10) 



where the cut correlator n(/i, Zcff (m)) reads 

ott(A') 



As is obvious from (jlOp . the cut correlator satisfies the equation 

dlogn(/^, Zcs{^J-)) 



d(l^) 



Eo. (12) 



The cut correlator n(/i, Zcff (/i)) is the quantity that actually governs the ex- 
traction of the ground-state parameters. 

The "fiducial" © range of /i is defined as the range where, on the one 
hand, the OPE reproduces the exact expression with better than some chosen 
accuracy (for instance, within, say, 0.5%) and, on the other hand, the ground 
state is expected to give a sizable contribution to the correlator. If we include 
only the first three power corrections (that is. Hi, 112, and Ha) we must require 
Lo/fi < 1.2. Since we know the ground-state parameters, we fix uj/n > 0.7, 
where the ground state contributes more than 60% of the full correlator. So 
our fiducial range is 0.7 < iv/^ < 1.2. 

We shall be interested in situations where the hadronic continuum is not 
known — which is typical for heavy-hadron physics and in the discussion of the 
properties of exotic hadrons. Can we extract the ground-state parameters? 

We denote the values of the ground-state parameters extracted from the 
sum rule pil| by E and R. The notations Eq and Rq are reserved for the exact 
values. In many interesting cases the ground-state energy may be determined, 
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Figure 1: The effective continuum threshold Zcs{fJ,) obtained by solving ([8]) for 
E = Eq and R = 0.7 Rq [long-dashed (blue) line], R = Rq [solid (red) line] and 
R — 1.15i?o [dash-dotted (green) hne]. 



e.g., from experiment. However, setting E — Eq does not help: still, for 
any R within a broad range, one finds a function Zce{fi, R) [Fig.[T] which solves 
the sum rule (jlOp exactly. Therefore, we conclude that in a limited range of 
fj, the OPE alone cannot say much about the ground-state parameters. What 
really matters is the continuum contribution, or, equivalently, Zcs{fi). Only by 
making some assumptions about Zoff (m) one is able to extract R. 

Typically, one assumes Zcsip-) to be constant and imposes some criteria 
to fix its value. Rigorously speaking, a constant effective continuum threshold 
Zeff (a*) = Zc = const is incompatible with the sum rule ^TU\i . Nevertheless, such 
an Ansatz may work well, especially in our HO model: As seen from Fig. [U 
the exact Zeff(/i) is almost flat in the fiducial interval. Therefore, the HO model 
represents a very favorable situation for applying the QCD sum-rule machinery. 

Now, how to determine Zc? A widely used procedure is to calculate 

p. X dlogH(^,Zc) 

i^(M,-c) = ^^^y^, (13) 

which now depends on due to approximating Zcff{n) by a constant. Then, one 
determines /xq and Zc as the solution to the system of equations 
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yielding z^ = 2.454(jj, ^io/uj — 1 [Fig. [2]. Finally, one takes the value R{iio, ^c) 
as the sum-rule estimate for the quantity R. The error of R is usually obtained 
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Figure 2: Constant effective continuum thresliold Zc- E{fi) for three different 
values of Zc (a) and the corresponding (b). 

by looking at the range covered by R{fJ., Zc) when one allows for a variation of /i 
within the fiducial range. Following this procedure, one obtains in our case a 
good central-value estimate: R/Rq = 0.96. Since R{fi, Zc) is extremely stable in 
the fiducial range, one expects its true value to be rather close to the extracted 
value and, accordingly, assigns a very small error to the sum-rule estimate. 

Note, however, a dangerous point: (i) a perfect description of n(/i) with an 
accuracy better than 1%, (ii) a deviation of Zc) from i?o at the level of only 
1%, and (iii) an extreme stability of in the entire fiducial range conspire to 
lead to a 4% error in the extracted value of R\ Clearly, this error could not be 
guessed on the basis of the other numbers obtained, and it would be wrong to 
try to estimate the error from, e.g., the range covered by R when varying the 
Borel parameter /i within the fiducial interval. 



Let us summarize the lessons we have learnt from the above investigation: 

1. The knowledge of the correlator to any accuracy within a limited range 
of the Borel parameter /i is not sufficient for an extraction of the ground-state 
parameters since rather different models for the correlator, generically of the 
form of a ground state plus an effective continuum, lead to the same correlator. 

2. Modelling the hadron continuum by a constant effective continuum threshold 
Zc allows one to determine the value of Zc by, e.g., requiring the average energy 

to be close to Eq in the region of stability of the sum rule. In the model 
under discussion this leads to a good estimate, R/Ro — 0.96, with almost 
/x-independent R. The unpleasant feature of this extraction procedure is that 
the deviation of R from Rq is much larger than the variations of E{fi) and 

over the fiducial interval of fi. In particular, it would be wrong to assign 
the systematic error on the basis of the range covered by R{p) when fj, is 
varied within the fiducial interval. This means that the standard procedures 
adopted in QCD sum rules do not allow one to control the systematic errors. 
Consequently, no rigorous systematic errors for hadronic parameters extracted 
by sum rules can be provided. Let us also stress that the independence of 
the extracted values of the hadron parameters from the Borel mass ^ does not 
guarantee the extraction of their true values. 

Finally, in the model under consideration sum rules provide a rather good 
estimate for Rq, even though its error cannot be determined on the basis of the 
standard procedures adopted in sum-rule analyses. This may be a consequence 
of the following features of the model: (i) a large gap between ground state and 
the first excitation contributing to the sum rule; (ii) an almost constant exact 
effective continuum threshold. Whether or not the same good accuracy may be 
achieved in QCD, where the features mentioned above are absent, is not obvious 
at all and requires more detailed investigations. 

We would like to point out that with respect to the problem of assigning 
systematic errors to the extracted hadron parameters, the method of QCD sum 
rules faces very similar problems as the application of approaches based on the 
constituent quark picture: for instance, the relativistic dispersion approach 
yields very successful predictions for the form factors of exclusive D decays and 
provides many predictions for the form factors of weak decays of B mesons . 
However, assigning rigorous errors to these predictions could not be done so far. 



Acknowledgements. Wc thank R. Bcrtlmann, B. Grinstcin, and B. Stcch for 
interesting discussions and inspiring comments, and the Austrian Science Fund 
(FWF) for financial support under project P17692. 



References 

1. M. Shifman, A. Vainshtein, and V. Zakharov, NucL Phys. B 147, 385 
(1979). 

2. W. Lucha, D. Mehkhov, and S. Simula, "Systematic errors of transition 
form factors extracted by means of light-cone sum rules" , larXiv:0712.0178[ 

3. W. Lucha, D. Melikhov, and S. Simula, Phys. Rev. D 76, 036002 (2007); 
Phys. Lett. B 657, 148 (2007). 

4. W. Lucha, D. Melikhov, and S. Simula, Phys. Rev. D 75, 096002 (2007); 
W. Lucha and D. Melikhov, Phys. Rev. D 73, 054009 (2006); Phys. Atom. 
Nucl. 70, 891 (2007). 

5. V. Novikov, M. Shifman, A. Vainshtein, and V. Zakharov, Nucl. Phys. B 
237, 525 (1984). 

6. A. L Vainshtein, V. I. Zakharov, V. A. Novikov, and M. A. Shifman, Sov. J. 
Nucl. Phys. 32, 840 (1980). 

7. V. A. Novikov et ai, Phys. Rep. 41, 1 (1978); M. B. Voloshin, Nucl. Phys. 
B 154, 365 (1979); J. S. BeU and R. Bcrtlmann, Nucl. Phys. B 177, 218 
(1981); Nucl. Phys. B 187, 285 (1981); V. A. Novikov, M. A. Shifman, A. I. 
Vainshtein, and V. I. Zakharov, Nucl. Phys. B 191, 301 (1981). 

8. A. Le Yaouanc et ai, Phys. Rev. D 62, 074007 (2000); Phys. Lett. B 488, 
153 (2000); Phys. Lett. B 517, 135 (2001). 

9. D. Mehkhov and S. Simula, Phys. Rev. D 62, 074012 (2000). 

10. M. Jamin and B. Lange, Phys. Rev. D 65, 056005 (2002). 

11. D. Melikhov, Phys. Rev. D 53, 2460 (1996); Phys. Rev. D 56, 7089 (1997); 
Eur. Phys. J. direct C4, 2 (2002) | hep-ph/0110087 ; D. Melikhov and S. 
Simula, Eur. Phys. J. C 37, 437 (2004). 

12. D. Melikhov and B. Stech, Phys. Rev. D 62, 014006 (2000). 



